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ABSTRACT 


From presently available observations one can Infer that the Alf- 
v6nic turbulence measured In the solar wind, predominantly on trailing 
edges of high speed streams, is a mixture of modes with two different 
polarizations, namely, Alfv6nic modes and modes which are the incompres- 
sible limit of slow magnetosonic waves. Using Helios 2 magnetic data 
and a variance analysis, we have separated parallel (to the mean field) 
and perpendicular components of the fluctuations and studied the possi- 
ble correlation between such components which would be predicted as a 
consequence of the incompressible character of the turbulence. Correla- 
tions between eigenvalues of the variance matrix are also investigated 


and discussed. 
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INTRODUCTION 


In a recent critical analysis of the present observations of 
almost incompressible Alfvfinic fluctuations in the solar wind [Dobro- 
wolny et al . , 1980a], several noticeable properties of such fluctua- 
tions were pointed out. One of them, obtained by comparing observatio- 
nal results with the equations of incompressible magnetohydrodynamics, 
is that such fluctuations should be a mixture of modes with two diffe- 
rent polarizations, namely, purely Alfvfinic modes and modes which are 
the incompressible limit of slow magnetosonic waves. 

Furthermore, by regarding the observed fluctuations as a turbulen- 
ce (and not as simple waves), a model ">i limum variance matrix was con- 
structed reproducing some observed statistical properties of the fluc- 
tuations under quite general assumptions on the distribution of wave 
vectors in the turbulence. 

In this paper we use extensive sets of Helios 2 magnetic data, 
exhibiting 'pure' (i.e. almost incompressible) Alfvfinic fluctuations 
in the trailing edges of several streams to study directly the polari- 
zation state of such turbulence through correlations between components 
of the fluctuations parallel and perpendicular to the background magne- 
tic field. We also test same scaling relations between eigenvslues 
which are predicted by the turbulence model recalled above. 

We start recalling in section 2 some main points in the theoreti- 
cal analysis in [Dobrowolny et al., 1980a] on the nature of Alfvfinic 
fluctuations, which lead to the tests on the data which are later pre- 
sented. 

In section 3 we present some results of the minimum variance ana- 
lysis of Helios 2 magnetic data, referring to the trailing edges of 
several high speed streams. Such analysis evidences the presence of 
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Alfv6nic fluctuations with statistical properties similar to those 
already obtained in various previous investigations [Belcher and Da- 
vis, 1971; Burlaga and Turner, 1976]. 

In section 4 we study the correlation between the component <Bg 
of the fluctuations parallel to the average magnetic field and the 
larger perpendicular component <B^ . We recall that, in a previous work 
[Sari and Valley, 1976] several phase and amplitude correlations 

between Fourier components of 6B# and 6B X were presented (in, approxi- 

-4 -2 

mately, the frequency range 10 - 10 Hs). In this paper, the correla- 

tion is studied on the data themselves (at each instant of time within 
our resolution) and not between their Fourier components. 

In the same section 4 we also investigate the correlation between dif- 
ferent eigenvalues of the variance matrix. 

Section 5 summarizes and discusses the conclusions of this work. 

« 
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2. THEORETICAL CONSIDERATIONS ON THE NATURE 
OF INCOMPRESSIBLE MHD TURBULENCE 
IN THE SOLAR WIND 

The AlfvAnic turbulence associated with trailing edges of high 
speed streams in the solar wind [Coleman, 1968; Belcher and Davis, 

1971; Burlaga and Turner, 1976] can, to a good approximation, be re- 
garded as incompressible. This statement follows quite clearly from 

a comparison between the typical level of such fluctuations, <B / | <B> | ^ 

r 

0.310.4, with the estimated level of fluctuations in the magnitude of 

B, which is € I B ] / 1 <B> | *'•0 .03*C . 06 [Burlaga and Turner, 1976] and with 

the recent results on associated density fluctuations [Neugebauer et 

al., 1978], giving 6p/o''-$ I v | /V , where V is the AlfvSnic speed in the 

“A A 

average magnetic field <B> . 

It has been shown in [Dobrowolny et al., 1980a] that the observed corre- 

i 

lations between velocity and magnetic field fluctuations in regions of 
Alfvlnic turbulence imply that such turbulence is in a very remarkable 
state characterized by the absence of non linear interactions . Physical 
arguments, Dased on the structure of the MHD equations, have been 
further given [Dobrowolny et al., 1980b], indicating that the tendency 
toward a state without non linear interactions may be a general proper- 
ty of th ; development of incompressible MHD turbulence and not a pecu- 
liar one of the solar wind turbulence. 

If non linear interactions are absent, the fluctuations which 
are present can be expanded (whatever their amplitudes) as a linear 
Fourier series, i.e. 

„ i(K*r - w t) 

8i - 

«B(r,t) = — l «b(K) e (1) 

K 
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with w - K*V . On the other hand, the equation V*4B-0, which muat hold 
K ™ “A 

for each Fourier component, indicates that, for each K, there are two 
independent polarisation for the fluctuations 4b(K) (with 4b A K) , i.e. 


4b(K) i 4b^(K) 


+ 4b (K) s 
”2 


5 4b 1 (K) 


« 1 (K) + 4b 2 (K) I (K) 


(2) 


with the two independent directions and given by 
K x <B> 

I (K) * — 

|K x <B>| (3) 

V £1 ■ i§ x VP 

i 

Clearly, the coaponent 4b, (K) of the fluctuations parallel to e, is 
| -1 " / 1 

entireljy perpendicular to <B> and, hence, in terms of proper inodes of 

i ~ 

the incompressible plasma, is what is called properly an Alfv6n wave, 
whereas' the component 4t^(K), parallel to $ 2 * has also a non zero 
projection 6b# parallel to <B> and can therefore be said the incompres- 
sible limit of the slow magnetosonic mode [Ginzburg, 1970]. 

The existence of a component of fluctuations parallel to <B> is, 
on the other hand, clearly indicated by the statistical analyses of the 
observations [Belcher and Davis, 1971; Chang and Nishida, 1973; Burlaga 
and Turner, 1976; Bavassano et al. , 1978], through the method of mini- 
mum variance. Their results show that the minimum eigenvalues (X^) of 
the variance matrix (of eigenvalues X^X^X^), although smaller than 

the other two (X «(X ,X )), is however different from sero. 

3 12 
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A second consequence of the absence of non linear interactions, 
if one adds to it the hypothesis of incompressibility (p>const), has 
been shown to be [Dobrowolny et al. ( 1980a] the constancy of the magne- 
tic field magnitude 


B 


2 


const 


( 4 ) 


Eqs. (2) and (4) imply in turn a precise relation hscween parallel 

2 

and perpendicular components of the fluctuations, the power <B± in per- 
pendicular fluctuations having a parabolic dependence from 5B W . If we 
make the further hypothesis 


I «B^| /6Bj. < 1 


(5) 


which is consistent with the observational evidence of a minimum eigen- 
value (and will be seen to be most often verified on the data), we 
obtain from (4) 

6B /7 1 «B X 

TiT - T^FI* (6) 


with C a constant given by 


C 




i] 


(7) 


In (6), 6B„ refers obviously to the component at polarization e of the 

2 2 

fluctuations, whereas the 6B± on the right hand side includes both pola- 
rizations e„ and S . 

1 2 
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The presence of two polarizations in tha fluctuations* expressed 
by (2) and (3)* and tha constancy of magnetic field magnitude (4) have 
consequences on the properties of the variance matrix characterising 
such fluctuations. 

In Dobrowolny at al. [1980a] a model calculation of such varian- 
ce matrix waa performed for a turbulence characterised by (2) and (4). 
For the simplest possible model of turbulence, i.e. a statistically 
homogeneous and normal turbulence (triple correlations all equal to 
zero), and using the further hypothesis 


«b./«b < 1 

4- 1 


( 8 ) 


the results obtained give the following simple scaling relations bet- 
ween eigenvalues 

(»> 

* ‘l 

the first of which can also be understood, in dimensional terms, from 
eq. (6), (which in fact indicates that the power in parallel fluctua- 
tions canj be expected to be proportional to the square of the power 
in perpendicular fluctuations. 
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3. MAGNETIC FIELD TATA ANALYSIS 


The data analyaia performed ha; used the magnetic field measu- 
rements of Helios 2. This spacecraft, launched on January 15, 1976, 
has been injected in a solar orbit having an aphelion of 0.98 AU and 
a perihelion of 0.29 AU, with an orbital period of about six months. 
During the primary mission of Helios 2 (Janusry to April, 1S76) we 
have selected, by visual inspection of hourly averages of the solar 
wind data of the Max Plank plasma experiment on Helios 2 (given to 
us by H. Rosenbauer and h. Schwenn; see also Schwenn et al. [1977]), 
a high velocity stream which is observed by the spacecraft during three 
successive solar rotations at different distances from the sun. The 
three observations of the stream begin on days 48 (February 17), 74 
(March 14), and 103 (April 12) of 1976, at heliocentric distances of 
about 0.89, 0.68 and 0.3}. AU respectively. 

In the trailing edge of this stream we have performed, using 6 s 
averages of the magnetic field vector, a sistematic analysis of the 
eigenvalues and eigenvectors properties of the variance matrix of the 
magnetic field evaluated over time intervals of different duration, 
ranging from about 3 minutes to 3 hours. 

In Figure 1 we present distributions of the ratios X /X , X / X , 
2, 2 1 3 | 

°B '*i # *1**2**3 bcing eigenvalues of the variance matrix and 
the variance of the field magnitude. Also shown is the histogram of the 
angle between minimum variance direction and average magnetic field. 
These results have been obtained for the two day period 50-51 within 
the trailing edge of the stream and using 22,5 minutes as a time basis 
for the statistics. 
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Theae results, end others which we do not report here, concerning 
eigenvalue distributions and distributions of th* minimum variance 
direction over the other streams and with different time basis, give 
the typical statistical properties already obtain K y ',?veral authors 
[Belcher and Davis, 1971; Chang, 1973; Burlaga find Turn* o, 1976 Cavas- 
sano et al., 1976 j for almost incompressible Aifvdnic turbulence in 
the solar wind. 

A discussion of such statistical properties, in terms of the diffe- 
rent frequency ranges for the fluctuations and as a function of distan- 
ce from the sun, will be the object of a forthcoming paper [Bavassa- 
no et el., I960]. 

The existence of correlations between parallel and perpendicular 
components of the fluctuations as well as of possible correlations 
between eigenvalues has been investigated for those periods, of each 
trailing edge, for which the Alfvdnic characteristics were more evi- 
dent. To have homogeneous sets of data, che duration of these periods 
was furthermore chosen so as to account for the rotational velocity 
of the spn as seen from Helios at the different times of the observa- 
tions (as a consequence the periods chosen range from about 2 days 
to little less than 4 days). The periods investigated have been divided 
in intervals of 22.5 minutes and for each interval the average field 
has been evaluated, then the fluctuations of the individual measure- 
ments (6 s averages) around the mean have been resolved in terms of 
parallel and perpendicular components with respect to the average field 
direction. Eigenvalues of the variance matrix have also been computed 
for each interval. 
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4, RESULTS 


A first aim of the present work is to present experimental eviden- 
ce of the relation (6) between parallel and perpendicular fluctuations. 
This evidence does, on the one hand, prove the gross validity of condi- 
tion (4) of constancy of the magnetic field and, on the other hand, 
indicate that the incompressible MHD turbulence observed is a mixture 
of modes with an Alfv6nic and a slow magnetosonic wave polarization 
as indicated in (2) and (3). 

To have a first insight about the relation between parallo 1 and 
perpendicular fluctuations we may examine scatter plots of the varia- 
bles 


x 


6B^ 

I <B> | 2 


y * | <B> | 


( 10 ) 


Some examples of them, for different intervals of 22.5 minutes, 

are given in Figure 2. If relation (6) is verified we should obtain 

points distributed along a straight line with a slope of -0.5. As it 

is seen the experimental data tend indeed to show such behaviour. The 

values of the correlation coefficient, given in the figure, strongly 

indicate the existence of a linear relation between x and y with an 

angular coefficient close to the value expected from (6). 

The examples in Figure 2 are not particular ones but represent 

inste&u a quite typical situation. This is proved from Figure 3 where 

we give frequency distributions for the correlation coefficient R and 

for the quantities C and C , which are the angular coefficient of 

1 o 

the best fit straight line and its intercept with the y axis, respecti- 
vely. The quantities R, C and C are evaluated on the basis of 22.5 

1 o 
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minutes intervals. For each parameter we show three different histo- 
grams, corresponding to the three selected periods of observation of 
the high velocity stream at differeht%hei T iocentric distances. From 
the left and eervkral pahels it appears clearly that in the great majo- 
rity of cases: a) the -correlation coefficient is close to -1, indica- 
ting in fact #%trong consistent tendency to a linear relation between 
x and y; b) the slope of the best fit straight line has a broad peak 
around -0.5, which is the value expected. These two properties remain 
pfhct*f tally unchanged when the distance from the sue varies. On the 
contrary the hiptoggiffe in the right panel reveal a general increase 

in the values of C approaching the sun. On accoQnt of (7), this can 
e ■ 

be takeF as ah indication of an increasing relative amplitude of the 
fluctuations when the skin dis&E&nce decreases. 

To gain sonte indications about; the properties of ie individual 
fluctuations we sho^w in Figure 4 the overall frequency distributions 

for | / | <H> | and | 6 | / 6 Bj_ . These histograms are obtai- 

ned from all the values of the parallel and perpendicular components of 
the fluctuations as evaluated, starting from 6 s averages, in the va- 
rious 22.5 minutes intervals of the different periods analyzed. It can 
be seen that approaching the sun the relative amplitude of the fluctua- 
tions Increases, as noted before, and that the parallel contribution to 
the fluctuating field becomes more important ' * should keep in mind, 
however, that we are neglecting possible temporal variations in the 
turbulence associated with the high speed stream. 

For a small but not negligeable portion of data points (see histo- 
gram of | 6 R /7 j / 6 in Figure 4)' the condition | 6 J / <5 < 1 is net well 

satisfied, which leads to a parabolic curve in the plane x-y. However, 
an examination of many examples of scatter plots shows that most of the 
experimental points fall in fact in a branch of the parabola where the 


linear approximation is good, so that the value of our linear come- 
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lation coefficient remains high. The few points for which the condition 
| 6B// ! / $Bj_< 1 is not valid increase just slightly the slope of the best 
fit straight line, as it is seen both in Figure 2(c) and in the histo- 
gram of (Figure 3). 

In conclusion we can say that in most of the analyzed time inter- 
vals there is a strong tendency to a linear relation between 6B/// 1 <B> | 

2 2 

and 6Bjl/|<B>| (correlation coefficient close to -1). Furthermore the 
slope of the best fit straight line has in general values near to that 
expected from relation (6). In other words the experimental data seem 
to prove the gross validity of the theoretical relationship. 

As a second task in the paper we have investigated correlations 
between eigenvalues of the variance matrix, to test the degree to which 
the scaling relations (9 ) , predicted by a model recalled in section 
2, are verified, 

2 

Figure 5 shows an example of scatter plots of X^-X^ and ^ -X^ 

(all eigenvalues being conveniently normalized), with superposed the 
best fit straight lines and the computed correlation coefficients. The 
Figure refers to the days 76 and 77, when the spacecraft is observing 
the trailing edge of the stream at a distance of ^0.65 AU from the sun. 

Similar plots have also been produced from the other data sets 
utilized so far (at different distances from the sun), taking overall 
periods of different length and, as already mentioned, with a time 
basis different from 22.5 minutes for the statistics. 

In general, the results obtained are very similar to the ones 
presented in Figure 5 and the values of the correlation coefficient 

reported there (R = 0.67 for the X -X correlation and R = 0.59 for the 

2 12 
X -X correlation) can be said to be typical. 

J 1 

Although these values for the correlation coefficient are signifi- 
cant, we see that the scaling relations (9) appear to be verified in 
a much poorer way than relation (6) between parallel and perpendicular 
components of the fluctuations at individual data points. 
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5. CONCLUSIONS 



From the results presented in section 4 we can draw the conclusion 
that the relation (d) between parallel and perpendicular components 
of the MHD fluctuations, observed in trailing edges of high speed 
streams, appear to be extremely well verified on an extensive set of 
data. We must stress that, although (6) is a consequence of the cons tan* 
cy of magnetic field magnitude (4), the test which has been done is 
more general than a test of constancy of |B|. As it should appear clear- 
ly from the discussion in section 2, it is also a test on the polariza- 
tion state of the observed MHD turbulence which, in this way, is in 
fact demonstrated to be a mixture of purely AlfvSnic (perpendicular) 
modes and modes which are the incompressible limit of slow magneto- 
sonic waves. The problem of how such a state is reached, starting, 
for example, from an initial production of both compressible and incom- 
pressible MHD fluctuations in the solar corona, is a different one 
and quite open at present. 

As far as the correlations between the eigenvalues of the variance 
matrix are concerned, they, in general, tend to reproduce the scalings 
(9), predicted by a model calculation. However, the correlation coeffi- 
cients obtained, although significant, are not as high as those of 
1 2 

the 6B// - 6 B a correlation between individual data points. As the sca- 
2 

ling X « A follows, in a dimensional way, from the same relation (6) 

2 

between and 6B^ which turns out so well tested with the data, we 
must ask ourselves which can be the reasons of the poorer correlation 
found between eigenvalues. There may be indeed two such reasons, both 
related to the special hypotheses of the model calculation in Dobro- 
wolny et al. [ 1980a]. 

The first one is the hypothesis (8) ab^/Sb^ < 1 used in the calcu- 
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lations. Although the parallel component 6B// of the fluctuations at 
polarization e^ is indicated by the analysis (see Figure 4) as smal- 
ler, in the majority of cases, than the perpendicular component fiB Xt the 
perpendicular fluctuations come from two contributions, at the polariza- 
tions $ and e^, and there is no way to separate these from the data, 

so that we do not know a typical value for fib /fib and the hypothesis 

2X 1 X 

used (fib^/fib^ < 1) may be wrong. On Lhe other hand, the results of 
Dobrowolny et al. ( 1980a] would certainly change if one was able to 
repeat an analogous calculation by maintaining fib / fib^ * 1. 

A second important point is that the calculation leading to (9) 
was assuming normal distributions for the turbulent fluctuations. Whe- 
reas this was done for the sake of having a tractable example, there 
is no reason to expect the real MHD turbulence to verify such an hypo- 
thesis and, in general, by ; analogy with developed hydrodynamic turbulen- 
ce, which is essentially not normal [Orszag, 1970], we would expect 
also in MHD turbulence triple correlations between fluctuating quanti- 
ties to be different from zero. Going to the scheme of the calcula- 
tions in Dobrowolny et al . [1980a], one becomes easily convinced that, 
in the presence of triple correlations, the scaling (9) between eigen- 
values wiJl no longer hold true but would be replaced by more complica- 
ted relations. 

Thus the fact that the testing of (9) with the experimental data, 
presented in this paper, has not given such high correlation coeffi- 
cients as the testing of the single point relation (6), may also re- 
flect the deviation from normality of the MHD turbulence in the solar 
wind. 
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FIGURE CAPTIONS 


Figure 1, Statistical distributions for the two day period 50-51 of: 

(a) the ratio of the eigenvalues X^/X^; the ratio X^/X^; ^ t * ie 

angle 6 between the minimum variance direction and the average field 

vector (the range of 6 being divided in equal increments of cosfl to cor- 

2 

rect for solid angle effect); (d) the ratio o /X of the variance of 

B 1 

the field magnitude to the largest eigenvalue. 

The analysis has been performed on a time basis of 22.5 minutes. 

2 2 

Figure 2. Examples of scatter plot of the variables x=6B±/|<B>| and 
y=6B*/|<B>| for different intervals of 22.5 minutes. In each plot we 
give the initial time of the interval, the best fit straight line and 
the value of the correlation coefficient R. 

i 

Figure 3. Statistical distributions of: (a) the correlation coeffi- 
cient R; (b) the slope of the best fit straight line; (c) the intercept 
of this line with the y axis. The three different histograms for each 
parameter refer to the three successive observations of the trailing 
edge of the stream. 

Figure 4. Statistical distributions of: (a) 6B X / | <B> | ; (b) I $B # | / 1 <§*> | ; 
(c) |6B^|/6B^ . The three histograms for each parameter refer to the 
three successive observations of the trailing edge of the stream. 

2 

Figure 5. Scatter plots of and A ^ “ A 3 ^ A 1 » A 2 an< * A 3 the 

variance eigenvalues X^ f X^ and X^ normalized to the average field 
magnitude) on a time basis of 22.5 minutes for the two day period 
76-77. In each plot the best fit straight line and the value of the 
correlation coefficient R are also given. 
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